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We study the entanglement dynamics of discrete time quantum walks acting on bounded finite 
sized graphs. We demonstrate that, depending on system parameters, the dynamics may be mono- 
tonic, oscillatory but highly regular, or quasi-periodic. While the dynamics of the system are not 
chaotic since the system comprises linear evolution, the dynamics often exhibit some features similar 
to chaos such as high sensitivity to the system's parameters, irregularity and infinite periodicity. Our 
observations are of interest for entanglement generation, which is one primary use for the quantum 
walk formalism. Furthermore, we show that the systems we model can easily be mapped to optical 
beamsplitter networks, rendering experimental observation of quasi-periodic dynamics within reach. 



I. INTRODUCTION 

Quantum walks (QWs) [T]-[5] have emerged as an in- 
teresting sub-field of quantum information processing [6 . 
Most existing work into QWs has focused on their appli- 
cability to implementing quantum algorithms [El HI 13 E] • 
However, several authors have also considered their abil- 
ity to generate entangled states [9 12 , which may subse- 
quently act as a resource for other quantum information 
processing protocols. In this paper we expand on previ- 
ous studies into the entanglement dynamics of QWs. We 
show that in some circumstances regular oscillations or 
elegant monotonic behaviour emerges, whereas for some 
system parameters, highly irregular, quasi-periodic be- 
haviour emerges in the entanglement metric of the sys- 
tem. We focus on bounded linear graphs with reflect- 
ing boundary conditions and show that reflecting bound- 
aries are key in the emergence of quasi-periodic dynam- 
ics. We primarily focus on single walkers, but also pro- 
vide some discussion on the case of multiple walkers and 
post-selection. The framework we present could easily be 
applied to analysing arbitrary graphs with any number 
of walkers. We pay special attention to the optical imple- 
mentation of such QWs. 

Carneiro et al. [13 previously examined in detail en- 
tanglement dynamics in discrete QWs acting on regular 
graphs. In this work we build on this by focusing on the 
quasi-periodic nature of these dynamics. While highly ir- 
regular dynamics emerge in some circumstances, the dy- 
namics are strictly non-chaotic since the system has no 
non-linearities. Wojcik et al. p3] previously gave con- 
sideration of the quasi-periodic nature of QWs. Here we 
focus on entanglement dynamics with static system pa- 
rameters and discuss in what circumstances regular pe- 
riodic or irregular quasi-periodic behaviour emerge. We 
find that even with fixed, non-periodic parameters, i.e. 
no time-dependent or periodic choices in the coins, quasi- 



periodic or perfect periodic dynamics may emerge. 

Considering bounded graphs is physically motivated 
as numerous experimental groups are beginning to ex- 
perimentally demonstate QWs [T5]-[22], but clearly infi- 
nite graphs are not possible, so considering the effects 
of boundary conditions is relevant to such studies. Fur- 
thermore, we discuss the mapping of our QW formalism 
to the specific case of photonic implementation, which is 
relevant to present optical experiments [2QH22] . 



II. BACKGROUND 

A QW comprises a system of positions (vertices in a 
graph), each of which may be occupied by a walker (e.g. a 
particle). A QW is defined by couplings between position 
states on some graph G. A position state may literally 
refer to a physical position, such as a spatial mode in op- 
tics, or more abstractly to a level in a higher dimensional 
system. The walker is an excitation, such as a photon, 
which may occupy a position, and is a bipartite system 
consisting of a position, x, and a coin value, c, |x,c). 
The coin is a parameter that indicates the direction the 
walker will take through the graph at the next step. The 
archetypal QW consists of two operations, the coin oper- 
ator and the step operator. The coin operator coherently 
specifies the coin value, which in general may be a super- 
position, while the step operator updates the position of 
the walker based on the coin state. This is in contrast to 
a classical random walk where the coin operator proba- 
bilistically specifies the coin value. 

The standard illustrative example of a discrete-time 
QW is that of a single walker evolving along an infinitely 
long linear graph structure. In this case we can define the 
coin and step operators as, 

C\x,±l) = {\x,-l)±\x,l))/V2 

S\x,c) = |x + c, c), (1) 
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where C is the Hadamard coin in this example. After t 
time steps the evolution of the system is given by (SCy. 
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In general the coin operator may vary at each time step 
in which case the evolution is of the form Hi^i ^^i- 

A. Walker operator formalism 

The simple illustrative example presented above is ap- 
propriate when a single walker is present and the graph 
has infinite linear structure (see for example Ref. [4 ). 
However that formalism becomes unwieldy for multiple 
walkers acting on general graphs. Thus we turn to the 
walker operator formalism presented by Rohde et al. [23] . 
Here a walker creation operator, k;(x,c)''', acts on an 
empty system |0), to create a walker with given position 
and coin values, 

^(x,c)'^|0) = |x,c). (2) 

In the photonic case |0) corresponds to the vacuum state. 
The state of a single walker system on graph G can be 
written in the form 

a,,^(z,j)^|0). (3) 

Evidently, with the interpretation that a walker is a par- 
ticle occupying distinct physical positions [37 , the state 
is always of a biased W-type [24 , i.e. a biased superpo- 
sition of a single excitation across a number of positions. 
Thus, the class of entangled states that may be gener- 
ated is rather restrictive [38|, and does not include some 
states that are of practical interest such as GHZ [25 and 
cluster (or graph) states However, W-states are 

nonetheless of practical interest since they are more re- 
silient against loss than other entanglement classes, and 
have applications in quantum cryptography [28] , telepor- 
tation [29-31 and computation [6 . 

In analogy with the conventional description of a QW, 
we define the evolution of the system using coin and step 
operators. The coin operator fiips a quantum mechanical 
coin which determines the direction of the walker. The 
step operator then uses the coin value to propagate the 
walker. 

C : w{x,c)^ ^ ^ A^^^w{x,j)^ 

S: w{xjy ^ w{j,xy, (4) 

Here Ux is the neighbourhood of position x in G, which 
in general may be a directed graph [39 , and A^^^ are 
local operators acting on the respective neighbourhood. 
Thus, the walker may only ever evolve to positions within 
its neighbourhood in the graph. An example of a walker 
graph is shown in Fig. [l] Importantly, the above for- 
malism allows for the modelling of irregular graphs, i.e. 
graphs where the degree varies at different vertices. Note 
the subtle difference in the definition of the step operator 
between the formalisms presented in Eqs. [l]and|4| In the 
former the step operator does not alter the coin state. 



whereas in the latter both the position and coin state are 
updated. This slight change in formalism is adopted to 
allow for arbitrary graphs to be described. 



^(2) ^(4) 




FIG. 1: An example walker graph with five position states. At 
each position i, a local operation A^^^ is applied to x and its 
neighbours. 

Once the coin operator is applied, the coin parame- 
ter can be interpreted as the next position of the walker, 
while after the step operation it can be regarded as a 
memory of the walker's previous position value. These 
properties are necessary to ensure unitarity of the oper- 
ations. In practise it is often somewhat artificial to con- 
sider the coin and step as separate operators; they can 
be merged into a combined evolution operator, 

E = SC : w{x, ^ (5) 

From this formalism one can express the allowed walker 
operator transitions as a graph. An example of the tran- 
sition graph for a two-position system is shown in Fig. [2j 

4(2) 

Q «;t(l,l) wH2,2) Q^'^ 




FIG. 2: (Colour online) The allowed transitions between 
walker operators for two position states. Vertices in red are as- 
sociated with position 1, while vertices in blue are associated 
with position 2. 



B. Global system evolution 

Based on the local neighbourhood operations we can 
construct the unitary operator acting on the entire sys- 
tem. For a single walker QW we can write a system of 
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coupled, discrete time-series equations as 

k 

This can be rewritten in matrix form as 
d{t^l) Ud{t), 

where 



(6) 



(7) 



(8) 



As an example, for a single walker, two-position system 
with balanced Hadamard coins. 
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the global evolution is given by 
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and for a single walker, three-position system by. 
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C. Optical implementation 



We now turn our attention to the specific case of op- 
tical implementation using single photons. In the optical 
case the described formalism is a little detached from the 
physical architecture as photons may have a well defined 
position, via distinct spatial modes, but have no intrinsic 
coin (i.e. memory) state To overcome this limitation 
we can map the walker operator formalism to a photon 
creation operator formalism by defining a function which 
maps the tuple {x, c} to a single number n which cor- 
responds to the photon's spatial mode, denoted by the 
photon creation operator aj^, where 



{x — l)d + c. 



(12) 



The reverse mapping can easily be performed using mod- 
ulo arithmetic. Having performed this mapping, it is 
straightforward to physically implement the walker oper- 
ator formalism in an optical context. Examples for map- 
ping two- and three-position walks are shown in Fig. |3] 
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FIG. 3: Mapping a two- (top) and a three-position walk (bot- 
tom) to a coinless photonic implementation. The coin opera- 
tor (C) applies the unitary evolution A^^"^ to each 'bundle' of 
modes corresponding to a given position i. The step operator 
(S) then simply permutes spatial modes appropriately. 



III. ENTANGLEMENT DYNAMICS 

Next we consider the entanglement dynamics of a QW. 
For a single walker we characterise the entanglement of 
the system using the Shannon entropy of the probability 
distribution of the walker across the position/coin basis 
states, 



E ■ 



■El 

ijeG 



'log2 



(13) 



This metric varies between 0, for a completely localised 
(i.e. unentangled) state, and log2 |Gp for a maximally 
entangled state (i.e. balanced W-state). 

For multiple walkers we can introduce the Meyer- 
Wallach entanglement measure [32 , given by 



E = 1 



1 



i=l 



(14) 
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where pi is the reduced state of subsystem i. £^ = for 
a completely separable system, and £^ > is a signature 
of entanglement. This measure effectively averages over 
the purity of the reduced state of each site in the system. 
This was previously applied to Q Ws in Ref. [12^ . 



A. Linear walks with a single walker 

We focus on the case of a single walker evolving 
through a bounded linear graph. We will study biased 
Hadamard coins, 



2. Biased Hadamard coins 

Next we turn our attention to the effects of coin bias 
in the evolution of the walk. Fig. [5] shows the entangle- 
ment and amplitude dynamics for a |G| = 5 linear graph, 
but with different coin biases. For (5 = the coin does 
not create a superposition of coin values, thus the walker 
wanders around while remaining completely localised at 
all times, and E = {). For 5 — 0.01 the behaviour becomes 
quasi-periodic but regular, and as (5 is increased further 
the dynamics become increasingly irregular over the il- 
lustrated time scale. 



(15) 



and examine the dependence of the entanglement dynam- 
ics on the coins' bias for varying graph sizes. For 5 = 1 
this operator reduces to the Pauli Z operator, and for 
6 = to the Pauli X operator, neither of which create 
superpositions. However for < (5 < 1 this coin creates 
superpositions of coin states, which subsequently creates 
entanglement. 



IV. PROPERTIES OF THE ENTANGLEMENT 
DYNAMICS 

We now examine some of the characteristics of the en- 
tanglement dynamics of the system. Some of the features 
are similar to chaotic systems, however our system is not 
truly chaotic since the evolution of the system is lin- 
ear. We consider two aspects of the system's evolution: 
(1) sensitivity to the system's parameters, and (2) the 
periodicity /quasi-periodicity of the dynamics. 



A. Sensitivity to system parameters 



1. Balanced Hadamard coins 

First consider the case of balanced Hadamard coins 
(i.e. 6 = 1/2) for linear graphs of varying lengths. Sev- 
eral examples of the entanglement generated on graphs 
of varying degrees are illustrated in Fig. |4] In the case 
of \G\ = 2, there are four basis states. The walker is ini- 
tialised in a localised state and E = 0. After one time step 
the state evolves into a Bell state across two of the basis 
states, after two time steps to a W-state across all four 
basis states (i.e. the maximum possible entanglement), 
after three time steps back to a Bell state, and then finally 
back to a localised state. Thus, the entanglement exhibits 
perfect periodicity with T = 4. For |G| = 3 the entangle- 
ment metric becomes quasi-periodic, with apparent irreg- 
ularity within the periods, and for \G\ =5 the behaviour 
seems to be completely irregular. For a very long graph, 
\G\ = 100, the entanglement initially increases monotoni- 
cally to an asymptote since the walkers are spreading out 
and able to occupy an increasing number of positions. It 
is not until well after the walker hits the boundary of 
the graph that irregular behaviour emerges. Again, af- 
ter this point the behaviour is quasi-periodic with highly 
irregular oscillations within the period. We note that in 
the case of \G\ = 100, the entanglement dynamics prior 
to reaching the graph's boundary mimic those observed 
in Ref. |9], where an infinite linear graph was considered. 
Clearly this is expected since prior to the walker hitting 
the boundary, the length of the graph is irrelevant to the 
dynamics of the system. 



There are several parameters characterising our QW 
system: (1) the graph structure, (2) the initial condition, 
and (3) the coin bias. In this case we will focus on a linear 
graph and consider the effects of coin bias. In Fig. [6] we 
consider a \G\ =5 linear walk with two slightly different 
coin biases, S = 0.5 and S = 0.51. Both walks have the 
same initial condition where the walker is initially com- 
pletely localised at the middle of the graph. Evidently 
the dynamics exhibit high sensitivity to the different coin 
bias parameters. Both evolutions initially exhibit similar 
irregular dynamics but quickly diverge and follow dif- 
ferent trajectories. We conclude that the entanglement 
dynamics in QWs are in general highly sensitive to the 
system's parameters. 

B. Periodicity 

Next we consider the periodicity of the system's evo- 
lution for different graph sizes. Revivals in QWs have 
previously been studied by Stefahak et al. [33 . Visually 
(see Fig.[5|, it is evident that the dynamics of the system 
are quasi-periodic, with repetitive structure over vary- 
ing time scales. To analyse this more quantitatively we 
can consider the eigenvalue spectra of the unitary matri- 
ces describing the evolution of the systems. Because we 
are dealing with a closed system (i.e. the magnitude of 
the wave- function is always preserved — it is never dis- 
sipative), the eigenvalues are all of unit magnitude, but 
in general complex. Furthermore, in general some of the 
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(d) \G\ = 100 
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FIG. 4: (Colour online) Entanglement and position amplitude evolution for linear graphs of various degree, 1^1 , and a fixed 
balanced Hadamard coin, 6=1/2. In the position amplitude plots, black corresponds to unit amplitude while white corresponds 
to zero amplitude. Blue lines show the entanglement dynamics, while red lines show the maximum possible entanglement (i.e. 
that of a balanced W-state). In (d) the insets show the monotonic and irregular/quasi-periodic sections in the evolution. Note 
the different time scales on the graphs. 



eigenvalues will be degenerate. The arguments of the dif- ferent eigenvectors of the unitary matrix. Thus the eigen- 
ferent eigenvalues then quantify the periodicity of the dif- value spectrum quantifies the periodicity of the system's 
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FIG. 5: (Colour online) Entanglement and position amplitude evolution for different coin biases, 6, and fixed graph length 
1^1 = 5. Note the different time scales on the graphs. 




FIG. 6: (Colour online) Comparison of the entanglement dy- 
namics for a 1^1 = 5 linear walk with 6 — 0.5 and S — 0.51, 
showing the system's sensitivity to slight perturbations in the 
coin bias. 



are rational, perfect periodicity will be observed, whereas 
when some of the arguments are irrational the system can 
be made to revive to within arbitrary accuracy over a suf- 
ficient timescale, although perfect periodicity will not be 
observed. Specifically, if the arguments of the eigenvalues 
can be expressed in rational form, Xj = e^^*^^ /^^ , then the 
periodicity of the system for a completely general state 



IS 



T = lcm(ti, . . . ,t|c'|2 



(16) 



where tj = qj / gcd{pj , qj) [41]. If we relax the conditions 
such that the system must revive up to a global phase 
then the periodicity may be a fraction of T. Obviously in 
special cases, such as where the walker is initialised into 
a linear combination of some subset of the eigenstates of 
/7, the period may be less than T. 



evolution over different time scales and as a function of 
the initial state. 

The eigenvalue spectra for different linear graph sizes 
are illustrated in the complex plane in Fig. ^ Note 
the bunching of the eigenvalues as \G\ is increased. For 
\G\ = 2 and \G\ = 4 all eigenvalues have arguments of ra- 
tional multiples of tt, whereas for |G| = 3 and |G| = 5 
at least some of the eigenvalues' arguments are irra- 
tional multiples of tt. Thus, the \G\ = 2 and |G| = 4 lin- 
ear graphs will exhibit perfect periodicity, whereas the 
|G| = 3 and |G| = 5 graphs will have infinite periodicity. 
More generally, when the arguments of all the eigenvalues 



V. MULTIPLE WALKERS AND 
POST-SELECTION 

Having studied the situation of single walkers, which 
we have shown lead to quasi-periodic behaviour, we now 
turn our attention to the case of multiple walkers. Multi- 
walker systems are of interest, particularly in the opti- 
cal context, since this is where exponential complexity 
in the system emerges and the possibility of exponential 
speedup of algorithms presents itself [34|. The power of 
the presented walker operator formalism is its ability to 
easily model higher numbers of walkers in a QW sys- 
tem. We show that with multiple walkers, perhaps not 
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FIG. 7: (Colour online) Eigenvalue spectra for linear graphs 
of different lengths 1^1. Note the bunching of the eigenvalues 
as 1^1 is increased. 



surprisingly, chaotic behaviour is still not possible, even 
with the addition of post-selection. It is well known that 
post-selection can induce non-linear effects in optical sys- 
tems [42] [35] , which is the basis of linear optics quantum 
computing (LOQC) [36^, thus opening the possibility of 
chaotic behaviour, since non-linearities are a necessary 
condition for the emergence of chaos. We explicitly work 
through the case of two walkers, but the expressions can 
easily be generalised to more walkers. 

For simplicity we will analyse the specific situation of 
photonic systems, since then we can ignore internal de- 
grees of freedom (i.e. coin states) and focus just on the 
position state, thereby simplifying the analysis. However, 
as discussed previously, a walker operator formalism can 
always be mapped to a photon creation operator formal- 
ism. Thus there is no loss of generality in restricting our- 
selves to the photonic case. 

Suppose we begin with two photons rather than one. 
Then the input state can be expressed as aja^, where i 
and j denote distinct optical modes. Then, after some 



unitary transformation of the creation operators, the 
state of the system will be of the form [23j 

k I 

kl 



(17) 



kl 



Thus the transformation is a bilinear map from the initial 
two photon state to a superposition of two photon states. 
By linearity, the coefficients of any superposition of differ- 
ent two-photon input states will be transformed linearly. 
Applying the same ideas as before, taking the eigenvalues 
of M will yield an eigenvalue spectrum where each of the 
eigenvalues gives a phase evolution of the correspond- 
ing eigenstate. Thus, quasi-periodicity may emerge. As 
an example comparison Fig. [Sj shows the Meyer- Wallach 
entanglement measure for the one- and two-walker cases 
with |G| = 5 and S = 1/2. Evidently, while distinct, the 
two graphs share many of the main features. Applying the 
same techniques as before, one could calculate the eigen- 
value spectrum of the unitary matrices to determine the 
overall frequency components in the dynamics. However, 
as the number of walkers increases the size of the unitary 
matrix grows exponentially as will the number of eigen- 
values. Thus, the dynamics will in general become more 
complex as the number of walkers increases. 

Next consider the situation where post-selection is per- 
formed. That is, we post-select on measuring a single 
photon at one output mode and consider the dynamics 
of the remaining photon. Upon measuring a single pho- 
ton at mode m the post-selected state of the system is 
given by [23] 

«]n {Uim Ujia\ + Ujm ^i/c^) = ^ijk^^l' 
I k k 

(18) 

Thus, we are left with a product state where one 
photon is projected into mode m and the remain- 
ing photon undergoes a linear map characterised by 

^ijk ~ ^imUjk + UjmUik' 



VI. CONCLUSION & FUTURE WORK 

In summary, we have considered the dynamics of an 
entanglement metric in discrete QWs. While we have fo- 
cused on bounded linear QWs with a single walker, the 
walker operator formalism adopted from Ref. [23 can 
easily be employed to generalise these results to arbi- 
trary graph structures with any number of walkers. Our 
results are particularly relevant to present day experi- 
ments, where linear walks are considered, but obviously 
infinite graphs are not possible. 

Our simple examples illustrate that highly complex, 
quasi-periodic dynamics can emerge even in small QW 
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|G| = 5 




boundary conditions are pivotal in the appearance of ir- 
regular dynamics, as similar studies on infinitely long 
lines demonstrate much more regular dynamics. Further- 
more, even with the introduction of multiple walkers and 
post-selection, non-linear evolution cannot emerge. 

We demonstrated that the walker operator formalism 
upon which our study is based, can easily be mapped to 
a linear optics photonic system. Thus, while challenging, 
experimental observation of irregular, quasi-periodic dy- 
namics in linear optics systems is realistic with present- 
day technology. Indeed, with the addition of bound- 
ary conditions, recently demonstrated experimental QWs 
[2QH22] could demonstrate such dynamics. 

In some previous studies it was observed that en- 
tanglement increases monotonically with time. However, 
our results indicate that in the bounded case this is not 
true and complex dynamics emerge. Thus, if the goal is 
to generate the maximal degree of entanglement, e.g. as 
a resource for other protocols, simply waiting for a very 
long time is not necessarily the best option. Therefore, 
timing is critical and understanding the precise nature of 
the entanglement dynamics is of importance. 

Future work should give consideration to potential uses 
for the generalised W-type entanglement that is produced 
by QW systems. Considering more elaborate graph struc- 
tures and alternate coin types would be valuable. This 
is something that can easily be approached using the 
walker operator formalism. In particular, the question of 
whether this can be used for efficient universal quantum 
computation or other quantum information processing 
applications is of importance. 



FIG. 8: (Colour online) The Meyer- Wallach entanglement 
measure for the one- (top) and two-walker (bottom) cases 
for a 1^1 = 5 linear graph with 6 — 1/2. 
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